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I.  Introduction 

Light  ion  beam  inertial  confinement  fusion  experiments 

being  conducted  on  PBFA  II  at  Sandia  National  Laboratories 

currently  employ  barrel  geometry  applied-B  ion  diodes^.  In 

diodes  of  this  geometry,  ions  are  produced  on  the  inner  surface 

of  a  barrel  shaped  structure  and  are  focused  over  a  distance  of 

2 

less  than  a  meter  to  the  centerline  of  the  barrel  .  The  close 

proximity  between  the  target  and  the  accelerator  structure  could 

lead  to  significant  radioactive  contamination  of  the  accelerator. 

One  possible  scheme  for  avoiding  chis  problem  is  to  use 

extraction  geometry  ion  diodes  to  focus  the  ion  beams  into  z- 

3 

discharge  plasma  channels  .  The  z-discharges  serve  as  transport 

channels  for  the  ion  beams.  They  can  provide  complete  charge  and 

4  5 

current  neutralization  for  the  propagating  beams  '  ,  and  because 
of  the  azimuthal  magnetic  field  created  by  the  discharge  current, 
the’ beam  ions  can  be  confined  to  propagate  within  the  channel 
radius^.  Transport  channels  that  are  several  meters  in  length 
can  provide  sufficient  standoff  between  the  target  and  the 
accelerator  to  protect  the  accelerator  from  the  deleterious 
effects  of  the  target  explosion.  In  addition  the  length  of  the 
transport  channels  allows  the  ion  beam  to  be  spatially  bunched  by 
the  use  of  a  carefully  programmed  accelerator  voltage  waveform'7. 
This  bunching  is  crucial  for  producing  the  power  level  and  short 
duration  pulse  needed  to  drive  the  target. 

A  problem  with  channel  transport  of  intense  ion  beams  is 
that  losses  associated  with  collisions,  instabilities,  and  the 
hydrodynamic  response  of  the  channel  to  the  passage  of  the  beam 
can  become  very  severe  in  small  radius  channels.  However,  these 

g 

losses  can  be  decreased  by  increasing  the  channel  radius  .  For 
example,  for  parameters  of  interest  to  the  Advanced  Pulse 
Experiment  (APEX)  of  Sandia  National  Laboratories  (30  MeV  Li*^ 
ion  beams  carrying  MA  level  currents  with  pulse  durations  of  10  - 

Manuscript  approved  August  3!,  1988. 
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50  ns),  channels  with  radii  of  2  to  3  cm  are  required  to 
sufficiently  minimize  the  beam  losses.  The  problem  with  large 
radius  transport  is  it  results  in  a  mismatch  betveen  the  beam 
radius  and  the  target  radius  vhich  is  1  cm  or  less. 

A  solution  to  this  radius  mismatch  problem  is  to  place  some 
type  of  ion  beam  focusing  cell  at  the  exit  of  the  transport 
channel.  A  proposed  focusing  cell  vhich  is  currently  under 
theoretical  and  experimental  investigation  is  the  one  eighth 

9 

betatron  wavelength  final  focusing  concept  .  In  this  concept  the 
focusing  cell  consists  of  a  short,  high  current  z-discharge.  As 
an  ion  beam  enters  the  focusing  cell,  the  sudden  increase  in  the 
azimuthal  magnetic  field  strength  over  that  in  the  transport 
channel  causes  the  ion  beam  to  pinch  inward.  At  the  exit  of  the 
focusing  cell,  the  majority  of  the  beam  ions  have  negative  radial 
velocities  so  that  the  beam  comes  to  a  focus  at  a  short  distance 
past  the  exit  of  the  cell.  The  beam  radius  at  the  focal  plane  is 
determined  by  the  ratio  of  the  transport  channel  discharge 
current  to  the  final  focusing  cell  discharge  current. 

The  natural  tendency  of  the  one  eighth  betatron  wavelength 

focusing  cell^  is  to  produce  a  focused  ion  beam  number  densi  ty 

profile  that  varies  as  1/r.  However,  the  .Initial  phase  space 

distribution  of  the  ion  beam  can  be  chosen  so  that  the  1/r 

profile  includes  an  off-axis  peak.  Figure  1  shows  an  example  of 

a  typical  1/r-like  number  density  profile.  This  profile 

corresponds  to  a  beam  that  was  transported  through  a  3  cm  radius, 

25  kA  discharge  current  z-discharge  and  focused  with  a  175  kA  one 

eighth  betatron  wavelength  focusing  cell.  The  ion  beam  was 

injected  into  the  channel  in  such  a  way  that  its  initial  radial 

phase  space  distribution  filled  a  rectangular  region  in  the 

radial  phase  space  plane  defined  by  0  <  r  <  2.12  cm  and 

2  2  1/2 

0  <  Iv  /v,  I  <  0.0598,  where  v  =  (v  +  v  )  .  Beam  ions  had  no 

1  r  d  1  d  r  z 

angular  momentum  (v^  =  0).  Figure  2  shows  an  extreme  example  of 
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a  number  density  profile  with  an  off-axis  peak.  The  difference 
between  this  example  and  that  in  Figure  1  is  that  ions  were 
injected  into  the  transport  channel  with  initial  radii  and  pitch 
angles  satisfying  1.94  cm  <  r  <  2.12  cm  and 

0.0548  <  |v  /vb|  <  0.0598.  In  this  second  example,  the  area  of 

radial  phase  space  covered  by  the  ion  beam  at  the  point  of 

injection  into  the  transport  channel  was  much  smaller  than  that 

-4 

in  the  first  example  (9.0  x  10  cm-radians  compared  to 
1.3  x  10'1  cm-radians).  Thus,  to  deliver  the  same  beam  current 
in  the  second  example  as  in  the  first,  the  ion  source  must  be 
much  brighter.  One  purpose  of  this  paper  is  to  show  how  the 
brightness  of  the  ion  source  must  increase  if  the  focused  beam  is 
to  have  a  nigh  concentration  of  ions  in  an  off-axis  peak.  In  the 
next  section,  the  theory  of  one  eighth  betatron  wavelength 
focusing  is  described  qualitatively.  Section  III  will  use  that 
theory  to  describe  how  the  initial  beam  injection  conditions  must 
change  if  the  focused  beam  is  to  be  concentrated  in  an  annular 
region.  These  results  will  be  coupled  to  two  simple  models  of 
extraction  ion  diodes  to  determine  the  ion  source  brightness 
requirements.  These  requirements  will  indicate  the  degree  of 
difficulty  involved  with  modifying  the  focused  beam  number 
density  profile  solely  by  modifying  the  initial  ion  beam  phase 
space  distribution.  In  Section  IV,  the  efficiency  with  ..iich 
an  ion  beam  is  delivered  to  any  given  annular  region  is 
calculated.  Section  V  summarizes  the  results  of  this 
investigation. 


II.  The  Theory  of  One  Eighth  Betatron  Wavelength  Focusing 

The  mathematical  description  of  one  eighth  betatron  wavelength 
focusing  of  intense  ion  beams  is  based  on  a  Vlasov  equation  analysis 
of  charge  and  current  neutralized  propagation  of  zero  angular  momentum 
ions  in  the  z-discharge  plasma  of  both  the  transport  channel  and  the 
final  focusing  cell^®.  While  the  mathematical  details  of  that  theory 
will  not  be  repeated  here,  the  qualitative  features  of  the  theory  will 
be  presented  along  with  the  results  that  will  be  required  in  the  next 
section. 

Figure  3  illustrates  the  main  components  of  the  beam  production, 
transport  and  focusing  system.  The  ions  are  produced  and  accelerated 
in  the  extraction  ion  diode,  which  also  focuses  the  ions  onto  the 
entrance  of  the  transport  channel.  The  transport  channel  is  a  long  z- 
discharge  plasma.  The  density  of  the  plasma  in  the  discharge  is 
sufficiently  high  that  the  ion  beam  is  completely  charge  and  current 
neutralized,  while  it  is  not  so  high  as  to  cause  a  large  amount  of 
collisional  energy  loss  as  the  beam  propagates.  The  azimuthal 
magnetic  field  associated  with  the  discharge  current  confines  the  beam 
ions  to  travel  within  the  discharge  radius  as  long  as  they  enter  the 
channel  at  sufficiently  small  radii  and  with  sufficiently  small  pitch 
angles.  The  final  focusing  cell  is  a  short  z-discharge  with  a 
discharge  current  several  times  higher  than  that  in  the  transport 
channel.  The  sharp  increase  in  the  magnetic  field  strength 
experienced  by  the  beam  as  it  moves  out  of  the  transport  channel  and 
into  the  focusing  cell  causes  the  beam  envelope  to  pinch  inward.  The 
length  of  the  focusing  cell  is  chosen  so  that  at  its  exit,  the  beam  is 
midway  into  its  pinching  phase.  The  beam  will  then  continue  to  pinch 
inward  as  it  propagates  through  field  free  space.  This  ballistic 
propagation  distance,  which  is  comparable  to  the  length  of  the 
focusing  cell,  is  considered  to  be  important  for  isolating  the  target 
from  the  plasma  in  the  focusing  cell. 
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The  dense  plasma  of  the  transport  channel  and  of  the  final 

focusing  cell  can  completely  charge  and  current  neutralize  the  ion 

^  3 

beam  ’  .  The  ion  beam  can  thus  be  regarded  as  a  collection  of  ioas 
propagating  independently  of  one  another  under  the  influence  of  the 
azimuthal  magnetic  field  produced  by  the  discharge  current.  Under 
this  condition,  a  single  particle  analysis  is  helpful  for  describing 
most  of  the  essential  physics  of  transport  and  focusing.  A  beam  ion 
in  the  transport  channel  feels  a  restoring  force  in  the  radial 
direction  due  to  the  radial  component  of  the  vxB  force.  This  causes 
the  ion  to  execute  an  approximately  sinusoidal  trajectory  called  a 
betatron  orbit  as  it  travels  down  the  length  of  the  channel.  The 
wavelength  of  this  trajectory,  i.e.  the  axial  distance  between  turning 
points  on  the  same  side  of  the  axis,  is  determined  by  the  ion  energy, 
charge,  and  mass,  the  channel  radius,  and  the  discharge  current.  In 
what  is  known  as  the  paraxial  approximation,  all  ions  have  the  same 
betatron  wavelength  given  by 


where  all  quantities  are  measured  in  cgs-es  units.  For  example,  a 
+3 

30  MeV  Li+7  ion  in  3  cm  radius  channel  carrying  100  kA  of  discharge 
current  has  a  betatron  wavelength  of  111  cm.  The  paraxial 
approximation  is  valid  for  ions  with  small  pitch  angles;  it  allows  the 
axial  component  of  the  velocity  in  the  radial  equation  of  motion  to  be 
replaced  by  the  ion  speed  v^,  which  is  constant.  Assuming  that  the 
magnetic  field  is  distributed  according  to  BQ(r)  =  (2I(,/crc)x(r/rc) , 
the  trajectory  of  an  ion  is  determined  by 
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(2) 


c  cr 
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c  cr  r 
c  c 
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In  the  paraxial  approximation,  this  is  the  equation  describing  a 
simple  harmonic  oscillator.  In  the  radial  phase  space  plane,  the 
motion  of  an  ion  traces  out  an  ellipse  defined  by 


r 


2 

+ 


(3) 


where  r  is  the  turning  point  of  the  ion.  Eq.  (3)  shows  that  for  the 
turning  point  to  be  equal  to  or  less  than  the  channel  radius,  the 
discharge  current  must  be  sufficiently  strong.  The  discharge  current 
required  to  confine  a  beam  ion  that  enters  the  channel  at  radius  r£ 
with  radial  velocity  vrE  can  be  found  by  setting  the  turning  point  in 
Eq.  (3)  equal  to  the  channel  radius  and  solving  for  I  .  This  gives 


I 

c 


2 


mc  vb 
2q 


rE/vb) 


2,  2 

Vrc 


(4) 


By  setting  r^  and  v^E  to  the  largest  values  that  characterize  the  ion 
beam  radial  phase  space  distribution  at  the  channel  entrance  this 
result  can  be  interpreted  as  a  lower  bound  on  the  discharge  current 
required  to  confine  the  ion  beam.  Note  that  this  relationship  can 
also  be  obtained  from  the  equation  of  conservation  of  canonical  axial 
momentum  without  invoking  the  paraxial  approximation.  When  the  ion 


1 


•i 


1 


i 

i 


6 


reaches  the  final  focusing  cell,  its  trajectory  in  radial  phase  space 
changes  from  the  ellipse  given  by  Eq.  (3)  to  a  different  ellipse 
defined  by  an  equation  of  the  same  form  as  Eq.  (3),  but  with  the  final 
focusing  cell  discharge  current  used  to  determine  the  betatron 
wavelength.  This  new  ellipse  is  narrower  in  radial  extent  and  broader 
in  the  radial  velocity  direction. 

The  behavior  of  a  beam  propagating  in  the  transport  channel  is 
determined  by  the  collective  behavior  of  *"he  individual  beam  ions.  In 
the  paraxial  approximation,  each  ion  travels  on  an  ellipse  in  the 
radial  phase  space  plane,  and  the  axial  distance  required  to  complete 
an  elliptical  orbit,  i .e.  the  betatron  wavelength,  is  the  same  for  all 
ions.  Thus,  whatever  initial  radial  phase  space  distribution  of  ions 
is  injected  into  the  channel,  it  reappears  at  intervals  of  the 
betatron  wavelength.  Between  these  reappearances  the  radial  phase 
space  distribution  rotates  on  a  fixed  ellipse  in  the  radial  phase 
space  plane.  If  the  paraxial  approximation  is  not  made,  then  the 
betatron  wavelength  of  an  ion  depends  on  its  turning  point.  Thus,  if 
the  initial  radial  phase  space  distribution  is  such  that  there  is  a 
spread  in  turning  points,  and  this  is  the  case  for  all  but  a  highly 
contrived  distribution,  there  will  then  be  a  corresponding  spread  in 
the  betatron  wavelengths.  Because  of  this  spread  the  ion  beam  radial 
phase  spp-s  distribution  will  be  a  quasiperiodic  function  of  the 
propagation  distance,  rather  tnan  being  a  simple  periodic  function  as 
in  the  paraxial  approximation.  The  initial  distribution  will  still 
reappear  at  intervals  of  the  betatron  wavelength,  but  at  each 
reappearance,  it  will  be  slightly  more  deformed  than  it  was  at  the 
previous  appearance.  After  many  betatron  wavelengths  of  propagation, 
this  gradual  deformation  causes  the  ion  beam  radial  phase  space 
distribution  to  effectively  fill  in  the  elliptical  region  in  the 
radial  phase  space  plane  defined  by  the  smallest  and  the  largest 
turning  points  in  the  ion  beam  distribution.  For  example,  the 
rectangular  initial  phase  space  distribution  defined  by  <  r  <  rv 
Vrl  -  lvrl  -  vr2  s^own  *n  FiSure  4  fills  in  the  area  between  the 
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ellipses  (denoted  by  the  dashed  curves)  defined  by 
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after  many  betatron  wavelengths  of  propagation.  Phase  mixing  makes 
the  ion  beam  radial  phase  space  distribution  independent  of  the 
propagation  distance.  It  is  assumed  in  this  analysis  that  the  ion 
beam  is  completely  phase  mixed  by  the  time  it  reaches  the  final 
focusing  cell. 

When  the  ion  beam  enters  the  final  focusing  cell,  the  phase  space 

ellipses  defined  by  Eq.  (5)  undergo  a  deformation  in  response  to  the 

abrupt  change  in  the  magnetic  field  strength.  If  the  magnetic  field 

has  the  same  distribution  in  the  final  focusing  cell  as  that  which  has 

been  assumed  for  the  transport  channel,  i.e.  BQ(r)  *  (2Ip/crc)x(r/rc) , 

where  1^  is  the  final  focusing  cell  discharge  current,  then  this 

deformation  is  a  rigid  c1_ckwise  rotation.  It  is  convenient  in 

describing  this  deformation  to  use  a  slightly  different  phase  space 

plane.  Let  T  *  (Aa/2it)(v  /v,  ),  where  Aa  is  the  betatron  wavelength  of 
pro  p 

a  beam  ion  in  the  final  focusing  cell  (Eq.  (1)  using  1^  in  place  of 

Ic).  The  phase  space  trajectory  of  an  ion  passing  through  the 

focusing  cell  is  a  circle,  rather  than  an  ellipse,  in  this  new  r-T 

phase  space  plane.  In  the  r-F  plane  the  elliptical  annulus  occupied 

by  the  beam  undergoes  a  rigid  clockwise  rotation  such  that  the  angle 

that  its  major  axis  makes  with  respect  to  the  r-axis  is  -2itz/Aa 

0 

radians.  Thus,  after  a  distance  of  one  eighth  of  a  betatron 
wavelength,  the  ellipses  have  rotated  -45  .  At  this  point,  the  number 
of  ions  with  negative  radial  velocities  is  maximized.  The  beam  then 
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enters  field  free  space  and  propagates  ballistically  to  the  target. 
While  it  does  so,  the  ellipses  in  the  r-T  plane  continue  to  rotate  and 
undergo  an  area  preserving  deformation  marked  by  a  progressive 
increase  in  the  eccentricity.  At  the  focal  plane,  the  ellipses  have 
rotated  a  total  of  -90  from  where  they  started  at  the  entrance  to  the 
final  focusing  cell,  and  the  ion  beam  occupies  the  area  in  radial 
phase  space  defined  by 
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f'Vf 

'V 
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Js 

,Vb, 

UjtJ 

(6) 


where  A  «  I,/I  -  1.  In  the  usual  case,  r_ ,  v  -,  and  the  channel 
t  c  l  rz 

discharge  current  are  chosen  so  that  the  maximum  turning  point  is 
equal  to  the  channel  radius.  This  makes  the  fullest  use  of  the 
available  magnetic  field  strength.  Thus,  the  beam  is  focused  to  a 
spotsize  of  [2/(2+A)  ]  ^^r  . 


To  this  point,  the  theory  is  independent  of  the  exact  form  of  the 
initial  radial  phase  space  distribution.  After  transport,  the  beam  is 
confined  to  an  elliptical  annulus  in  the  radial  phase  space  plane  with 
the  inner  and  outer  radii  of  this  annular  region  corresponding  to  the 
smallest  and  the  largest  turning  points,  respectively,  that  will  occur 
in  that  distribution.  This  is  illustrated  in  Fig.  4.  However,  to 
calculate  number  density  profiles,  the  exact  form  of  the  initial 
radial  phase  space  distribution  must  be  specified.  For  the  purposes 
of  this  calculation,  the  initial  radial  phase  space  distributions  that 
have  been  considered  are  the  family  of  rectangles  defined  by 
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r^  <  r  <  lvr]_l  <  !vrl  <  | vr2  i  *  Within  this  rectangular  region  of 

radial  phase  space,  ions  are  distributed  in  such  a  manner  as  to 
produce  a  number  density  profile  that  is  flat  betveen  r^  and 
Because  of  the  cylindrical  geometry,  this  means  that  the  initial 
radial  phase  space  distribution  is  more  heavily  weighted  at  larger 
radii.  This  initial  distribution  is  convenient  for  calculations  and 
is  also  a  useful  model  for  the  possible  distributions  produced  by  an 
extraction  ion  diode.  After  phase  mixing,  the  phase  space 
distribution  of  the  beam  depends  only  on  the  distance  from  the  origin 
of  the  phase  space  plane,  apart  from  a  1/r  factor  that  reflects  the 
cylindrical  geometry.  At  the  focal  plane,  the  distribution  can  be 
expressed  as  a  function  of  p,  where  p  is  defined  by 


2 

p 


(2->  a)  r2 
2  r 


+ 


2(1  ♦  A)  _2 
(2  ♦  A) 


(7) 


It  is  also  convenient  to  express  the  combination  ( A,/2Jt)  v^/v^  as  a  new 

variable  y.  In  the  r-Y  phase  space  plane,  analogous  to  the  r-T  phase  • 

space  used  for  describing  the  focusing  cell,  the  phase  space 
trajectory  of  an  ion  in  the  transport  channel  is  a  circle  rather  than 

an  ellipse.  There  are  then  two  results,  depending  on  the  magnitude  of 

2  2  1/2  2  2  1/2 
(y2  +  r^)  relative  to  •  The  two  result  correspond  to 

the  orientation  of  the  rectangle  that  defines  the  initial  radial  phase 

2  2  2  2 

space  distribution.  For  y^  +  r^  <  ^  +  y^,  that  rectangle  is  oriented 
horizontally,  and  the  focused  beam  phase  space  distribution  is 

t 


l 


* 
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For  Yj  +  >  *2  +  Yi»  the  initial  Phase  space  distribution  rectangle  is 

oriented  vertically,  and  the  focused  beam  phase  space  distribution  is 
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where  in  both  Eqs.  (8)  ‘and  (9),  N  represents  the  average  density  of  ions 
in  the  region  of  the  radial  phase  space  plane  occupied  by  the  initial 

phase  space  distribution.  For  both  equations,  if  N  is  the  total  number 

0  2 

of  ions  present  at  the  source  plane,  then  N  is  given  by  N  =  N  Xa/2n(r„  - 

2  O  p  Z 

rl^Y2  “  Yl^’  ^he  num^er  density  profile  of  the  beam  at  the  focal  plane 

is  then  obtained  by  integrating  either  Eq.  (8)  or  (9)  over  the  radial 

velocity.  The  complete  expressions  for  the  number  density  profile  at  the 

focal  plane  are  given  in  the  Appendix.  These  results  can  be  used  to  show 

that  any  off-axis  peak  in  the  focused  beam  number  density  profile  occurs 

between  radii  that  correspond  to  the  intersection  of  the  ellipse  defined 

by  the  phase  mixed  radial  phase  space  distribution  and  the  v  axis,  i . e . 

the  range  of  radii  given  by  Eq.  (6)  with  v  =  0. 
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III.  Estimates  of  Required  Source  Brightness 


As  discussed  in  the  previous  Section,  the  number  density  profile  of 
the  focused  ion  beam  depends  on  the  initial  radial  phase  space 
distribution.  For  the  model  initial  distribution  chosen  for  this 
analysis,  in  which  the  ions  occupy  a  rectangular  region  in  the  radial 
phase  space  plane,  the  focused  beam  density  profile  depends  only  on  the 
relative  positions  of  the  four  corners  of  that  rectangular  region.  In 
order  to  produce  a  focused  ion  beam  number  density  profile  that  has  a 
strong  off-axis  peak,  the  initial  distribution  must  be  strongly 
concentrated  in  a  small  region  of  the  radial  phase  space  plane  away  from 
r  *  0.  The  "narrower"  the  initial  distribution  is,  the  stronger  will  be 
the  off-axis  peak,  where  "narrow"  refers  to  the  span  of  distances  from 
the  origin  of  the  radial  phase  space  plane.  For  the  rectangular  initial 
distribution  used  here,  the  relevant  dimension  is  the  length  of  the 
diagonal.  To  deliver  the  same  ion  beam  current  with  a  narrower  initial 
distribution  requires  that  the  ion  source  concentrate  an  equal  number  of 
ions  into  a  smaller  area  in  the  radial  phase  space  plane,  i.e.  the  source 
brightness  must  increase  as  the  area  of  the  radial  phase  space  plane 
spanned  by  the  initial  distribution  decreases.  To  estimate  the  required 
source  brightness,  it  will  be  required  that  the  elliptical  annulus  that 
defines  the  phase  mixed  and  focused  ion  beam  radial  phase  space 
distribution  intersect  the  v^.  axis  in  a  range  of  radii  that  corresponds 
to  a  target  annulus  defined  by  r^n  <  r  <  r  ,  as  illustrated  in  Figure 
5.  By  using  the  transport  and  focusing  theory  discussed  in  the  previous 
Section,  this  requirement  can  be  translated  into  a  calculation  of  the 
source  brightness.  This  requirement  is  a  minimum  condition  for 
concentrating  the  beam  onto  a  target  annulus  but  does  not  account  for  the 
ions  which  miss  the  target  annulus  at  the  focal  plane.  In  the  next 
Section,  these  ions  will  be  accounted  for  through  the  calculation  of  a 
delivery  efficiency. 

As  mentioned  previously,  the  off-axis  peak  in  the  focused  beam  number 
density  profile  occurs  within  the  range  of  radii  given  by  Eq.  (6)  with 


v  »  0.  The  outer  radius  of  this  range  can  be  chosen  by  using  the 
appropriate  ratio  of  discharge  currents.  The  inner  radius  of  this  range 
is  determined  by  r^  and  v  It  will  be  assumed  that  ^  and  are 
chosen  such  that 


r 


2 

2 


+ 


(10) 


The  required  ratio  of  the  transport  channel  discharge  current  to  the 
final  focusing  cell  discharge  current  is  then  determined  by 


rout 


i  .e. 
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1  . 


(ID 

(12) 


From  Eq.  (6),  and  using  Eq.  (11)  to  eliminate  the  factor  2/(2+A),  r^ 
v  j  must  satisfy 


and 


(13) 


Eq.  (13)  then  determines  the  injection  conditions  in  terms  of  the 
location  of  the  off-axis  peak  in  the  number  density  profile  of 
the  focused  ion  beam. 

To  obtain  requirements  on  the  ion  beam  source,  specific 

injection  schemes  must  be  considered.  Tvo  such  schemes  will  be 

considered  here.  One  scheme,  illustrated  by  Figure  6, 

corresponds  to  the  conventional  extraction  diode  focusing  scheme. 

The  extraction  diode  is  designed  to  focus  to  a  point;  the  spot 

size  is  limited  by  microdivergence.  This  scheme  will  be  called 

pinched  focusing.  The  second  scheme  is  illustrated  in  Figure  7. 

In  this  scheme  the  extraction  diode  is  designed  to  produce  a  beam 

that  is  hollow  at  the  entrance  to  the  transport  channel,  with  the 

width  of  the  annulus  determined  by  microdivergence.  This  scheme 

will  be  called  annular  focusing.  Rather  than  dealing  explicitly 

with  the  radial  veloc:  „ies,  the  pitch  angles  and  will  be 

used  where  or^  =>  v^/v^  and  *  vr2/vb*  ^or  ^oth  focusing 

schemes  it  will  be  assumed  that  the  pitch  angles  at  injection  are 

determined  primarily  by  the  injection  geometry.  This  assumption 

is  valid  as  long  as  the  microdivergence  is  sufficiently  small. 

The  anode  of  the  extraction  diode  will  be  modeled  as  an  annulus 

with  outer  radius  R  _  and  inner  radius  R.  .  The  focal  length  of 

out  xn 

the  diode  will  be  denoted  by  F.  The  maximum  and  minimum  pitch 
angles  at  injection  can  then  be  expressed  as 


I 


I 
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(H) 


These  considerations  allow  the  source  current  density  required  to 
produce  a  given  ion  beam  current  to  be  calculated  as  a  function 
of  the  diode  focal  length  and  microdivergence.  This  will  be  done 
first  for  the  pinched  focus,  then  for  the  annular  focus. 

Pinched  Focus 

In  this  injection  scheme,  r^  =  0  and  r^  is  determined  by 
microdivergence: 

r2  .  F  A0  (15) 

The  anode  surface  is  modeled  as  an  annulus;  the  inner  and  outer 
radii  determine  both  the  area  availab*^  for  production  of  the  ion 
beam  and  the  minimum  and  maximum  injection  angles.  The  anode 
must  produce  the  required  ion  beam  current.  This  requires  a 
certain  anode  current  density  JA  given  by 
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Tt  Tr2  -  R2  1 

v  out  in; 


(16) 


where  Ig  is  the  ion  beam  current  in  the  diode.  Of  course,  in 
reality  the  anode  current  density  would  have  some  spatial 
structure,  both  in  radius  and  in  azimuth,  but  its  average  over 
the  anode  area  must  equal  that  given  by  Eq.  (16). 

The  maximum  injection  angle  is  determined  by  the  requirement 
that  ions  are  confined  to  the  transport  channel.  From  Eqs.  (10) 
and  (15),  the  maximum  angle  can  be  calculated  as 


“2 


-  F2(A6) 


(17) 


The  minimum  injection  angle  is  determined  by  Eq.  (13)  with  r^ 
*  0.  This  gives 


*1 


2rt  rin 

X a  r  c 

3  out 


(18) 


Using  the  relation  between  the  anode  radii  and  the  injection 
angles  given  in  Eq.  (14)  yields  the  equation  for  the  anode 
current  density 
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i>|i|2  1 

it  l  2itJ  F 


F2  (l  -  f2)r2  -  F2(A0)2 


where  f  *  r.  /r  ..  The  anode  current  density  can  be  minimized 
in  out 

with  respect  to  the  focal  length  F  yielding  the  optimal  focal 
length 


ll  -  f‘ 

7 Tie 


At  the  optimal  focal  length,  the  anode  current  density  attains  a 
minimum  given  by 


A, min 


Cl  -  f2): 


These  results  are  plotted  in  Figs.  8,  9,  and  10.  In  these 

figures,  and  in  all  of  the  figures  to  follow,  I_  =  3.33  MA  (30 

MeV  *^Li),  rc  =  3  cm,  and  rQUt  =  0.8  cm.  Two  different  values  of 

the  final  focusing  current  are  considered,  1  MA  and  300  kA.  For 

the  high  current  case,  Eq.  (12)  requires  that  Io  =  36.9  kA.  For 

the  low  current  case  I  =  11.1  kA.  In  Fig.  8,  the  anode  current 

c 

density  is  plotted  as  a  function  of  the  focal  length  for  several 
different  values  of  the  microdivergence.  In  this  plot,  the  final 
focusing  current  was  1  MA  and  r . n  =  6  mm.  In  Figures  9  and  10, 
the  minimum  anode  current  density  is  plotted  as  a  function  of 
microdivergence  for  several  different  values  of  r^n  for  the  high 
current  case  and  the  low  current  case,  respectively. 


Figures  11  and  12  show  how  the  inner  and  outer  anode  radii 

must  vary  with  microdivergence  to  achieve  the  minimum  current 

density  of  Eq.  (21).  Figure  11  corresponds  to  the  high  current 

case,  and  Figure  12  to  the  low  current  case.  Also  shown  in  these 

plots  is  the  optimal  focal  length.  The  anode  radii  were 

calculated  from  the  relations  R.  *  Fa.  .  R  =  Fa  using  the 

in  mm  out  max 

optimal  focal  length. 

The  minimum  current  density  of  Eq.  (21)  depends  on  the 

2 

microdivergence  quadratically ,  so  by  dividing  J.  .  by  (A0)  ,  a 

Af  min 

quantity  that  depends  only  on  the  transport,  focusing,  and  target 
parameters  results.  Figure  13  shows  a  plot  of  this  quantity  as  a 
function  of  rin  for  the  high  current  and  low  current  cases. 
Multiplying  this  quantity  by  the  ion  energy  will  give  the  power 
brightness  required  by  the  target.  Figure  13  shows  that  as  the 
beam  profile  becomes  progressively  more  annular,  the  required 
power  brightness  increases  rapidly,  spanning  nearly  an  order  of 
magnitude  as  rin  increases  from  0  to  6  mm,  '■hen  increases  more 
rapidly  as  rin  continues  to  approach  rQUt- 

Annular  Focus 

Figure  7  illustrates  the  annular  focus  injection  scheme.  In 
this  scheme,  a^  =  0.  The  diode  microdivergence  determines  the 
spot  size  at  the  transport  channel  entrance  through  r^  - 
r^  =  FA0.  The  maximum  injection  angle  is  determined  by  the 
requirement  that  ions  be  confined  to  the  channel.  The  minimum 
injection  radius  is  determined  from  Eq.  (13)  with  v  ^  =  0. 
Denoting  this  result  as  r^  =  fr^,  the  maximum  injection  angle  is 
given  by 


where,  as  before,  f  =  r.  /r 

m  out 


The  anode  surface  is  modeled  in  this  scheme  as  a  disc  with 
radius  RQut  which  is  related  to  as  indicated  in  Eq.  (14). 
required  anode  current  density  ran  then  be  written  as 


As  in  the  previous  case,  this  may  be  minimized  with  respect  to 
the  focal  length  F.  The  optimal  focal  length  is 


F„„,  -  Is.  f  3  .  i  r  T~T  1 

0pt  A0  L"4  f  *  72  J1  +  8  f  J 


At  this  optimal  focal  length,  the  required  anode  current  density 


V 
+  8f 


Figures  14  and  15  show  JA  as  a  function  of  microdivergence  for  several 


different  values  of  r^n  for  the  high  current  and  low  current  cases, 
respectively.  Figures  16  and  17  show  the  anode  radius  required  to  produce 
the  given  ion  beam  current  with  the  minimum  anode  current  density.  Also 
plotted  are  the  optimal  focal  lengths. 

The  minimum  anode  current  density  is,  as  in  the  previous  case, 

2 

quadratic  in  the  microdivergence.  Figure  18  shows  how  m^n/(&0)  varies 
with  r^n  for  both  the  high  and  low  current  cases.  Comparing  this  with  the 
similar  result  for  the  pinched  focus  shows  that  the  annular  focus  places 
much  more  stringent  requirements  on  the  extraction  diode.  For  the  annular 
focusing  scheme,  the  required  power  brightness  increases  by  nearly  two 
orders  of  magnitudes  as  r .  n  ranges  from  0  to  6  mm. 
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IV.  Calculation  of  Delivery  Efficiency 

The  calculations  presented  in  Section  III  vere  based  on  the  requirement 
that  the  bean  be  injected  into  the  channel  in  such  a  way  that  an  off-axis 
peak  in  the  focused  bean  number  density  profile  occur  betveen  tvo  given 
radii  r^n  and  r  .  From  Figure  2,  one  can  see  that  by  choosing  r^  and  v  ^ 
close  enough  to  r^  and  v^*  a  highly  peaked  number  density  profile  can  be 
produced.  However,  not  all  of  the  available  ions  are  concentrated  into  any 
given  annulus  at  the  focal  plane.  Thus,  while  the  conditions  on  the  source 
brightness  found  in  the  preceding  Section  are  useful  for  determining  how 
difficult  it  is  to  concentrate  the  beam  into  an  off-axis  peak,  they  are 
under-estimates  because  only  a  fraction  of  the  beam  is  delivered  to  the 
region  r^n  <  r  <  rout  as  illustrated  in  Figure  5.  If  one  requires  that  a 

certain  ion  beam  current  be  delivered  to  the  annular  region,  then  the 
estimates  given  in  the  preceding  Section  must  be  increased  to  account  for 
the  fraction  of  the  ion  beam  that  does  not  appear  in  the  annular  region. 

In  this  Section,  the  efficiency  with  which  an  ion  beam  is  delivered  to  the 
annular  region  r^n  <  r  <  rQut  is  calculated,  and  these  results  are  applied 
to  the  estimates  of  the  required  source  brightness  from  Section  III. 

The  delivery  efficiency  h  will  be  defined  as  the  number  of  beam  ions 
delivered  to  the  annular  region  divided  by  the  total  number  of  beam  ions 
delivered  to  the  focal  plane.  Thus,  h  is  just  the  ratio  of  the  integrals 
of  the  number  density 


n 


xout 

n^(r)  rdr 

rin _ 

xfoc 

n^(r)  rdr 

.0 


(26) 
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where  nb(r)  is  the  number  density  of  the  ion  beam  at  the  focal  plane. 
In  the  integral  in  the  denominator  of  Eq.  (26)  is  the  radius  of 

the  focused  ion  beam.  It  has  been  left  arbitrary  to  account  for  the 
possibility  that  the  ion  beam  is  either  under  or  over  focused  with 
respect  to  r  .  Due  to  the  functional  form  of  n^(r),  the  integrals 
appearing  in  Eq.  (26)  must  be  calculated  numerically. 

An  upper  bound  on  the  delivery  efficiency  can  be  obtained  by 
considering  the  zero  emittance  limit.  Suppose  that  all  of  the  beam 
ions  are  injected  with  the  largest  radius  and  pitch  angle  consistent 
with  the  requirement  that  the  ions  be  confined  to  propagate  within  the 
channel  radius.  The  phase  mixed  radial  phase  space  distribution  of 
such  a  beam  is  given  by 


f(r,vr)  =.  S(r2  +  (\0/2n)2(vr/vb)2  -  r2)  (27) 

i.e.  in  the  transport  channel,  ions  are  restricted  to  move  on  the 
2  2  2 

circle  r  +  r  -  rch  in  the  r-r  plane.  At  the  focal  plane,  this 
distribution  is  modified  to 


£('*v  *  «'2  *  <y2,t)2(vv2  -  (rr-sK) 


(28) 


i.e.  in  the  r-r  plane,  ions  occupy  the  circle  of  radius  [2/(2+A)]rc- 

The  distribution  of  ions  on  this  circle  is  independent  of  the  angle. 

Thus,  the  number  of  ions  that  appear  in  any  given  annular  region 

r.  <  r  <  r  .is  determined  by  the  fraction  of  the  circle 
in  -  out 

r  +  r  *  [2/(2+A)]r2b  that  passes  through  the  annulus  as  illustrated 
in  Figure  19.  Using  r ^oc  as  shorthand  for  [2/(2  +  A))rc>  the  delivery 
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efficiency  for  the  zero  emittance  limit  is 


^EL 


2 

ji 


cos 


-1 


-1 

cos 


(29) 


For  the  case  of  r^oc  =  rout  considered  in  the  previous  Section,  the 
zero  emittance  limit  delivery  efficiency  is  simply 


^EL 
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cos 


rin 


out 


(30) 


For  example,  the  zero  emittance  limit  delivery  efficiency  for 

r.  *  0.6  cm  and  r  .  »  0.8  cm  is  462.  The  zero  emittance  limit 
in  out 

delivery  efficiency  will  be  used  to  compare  with  the  delivery 
efficiency  for  the  more  realistic  cases  of  pinched  and  annular 
focusing. 


The  delivery  efficiency  has  been  calculated  numerically  for  both 
the  pinched  focus  and  annular  focus  cases  analyzed  in  Section  III.  As 
in  Section  III  the  injection  of  the  beam  into  the  transport  channel 
was  assumed  to  take  place  at  the  optimal  focal  lengths  as  defined  by 
Eq.  (20)  for  the  pinched  focus  case  and  by  Eq.  (24)  for  the  annular 
focus  case.  Figure  20  shows  the  delivery  efficiency  for  the  pinched 
focus  case  as  a  function  of  Also  plotted  in  Fig.  20  is  the  zero 

emittance  limit  delivery  efficiency.  This  result  shows  that  for  small 
values  of  rin  the  efficiency  is  close  to  the  zero  emittai.^e  limit,  but 
as  r.  approaches  r  „ ,  it  drops  rapidly  below  r>  .  Figure  21  shows 
the  delivery  efficiency  for  the  annular  focus  case.  The  delivery 
efficiency  for  the  annular  focus  case  is  quit^  similar  to  that  for  the 
pinched  focus  case. 
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These  results  can  be  used  along  with  the  calculation  of  the 
required  source  brightness  from  Section  III  to  determine  the  source 
brightness  required  to  deliver  a  specified  amount  of  beam  current  to 
the  region  r^n  <  r  <  r  .  This  is  done  by  dividing  the  source 
brightness  plotted  in  either  Fig.  13  or  Fig.  18  by  the  delivery 
efficiency  for  the  appropriate  case.  Figures  22  and  23  show  the 
results  of  this  calculation  for  the  pinched  focus  case  and  for  the 
annular  focus  case,  respectively.  In  these  calculations  it  vas 
required  that  3.33  MA  of  30  MeV  Li+^  ions  produced  at  the  diode 
(before  time-of-flight  bunching)  be  delivered  to  the  specified 
annulus;  the  total  beam  current  is  greater  than  this  amount  by  the 
reciprocal  of  the  delivery  efficiency.  These  results  indicate  that 
the  requirement  that  the  specified  beam  current  be  delivered  to  a 
given  annulus  places  very  stringent  requirements  on  the  ion  source, 
particularly  for  narrow  annuli.  As  shown  in  Figs.  22  and  23,  the 
required  source  brightness  increases,  for  example,  by  about  two  orders 
of  magnitude  as  r^  increases  from  4  mm  to  7  mm  for  a  target  annulus 
with  an  outer  radius  of  8  mm. 


V.  Conclusions 


The  effect  of  the  initial  radial  phase  space  distribution  of  an 

ion  beam  on  its  radial  number  density  profile  after  being  transported 

and  focused  vith  a  one  eighth  betatron  vavelength  focusing  cell  has 

been  investigated.  The  theory  of  one  eighth  betatron  vavelength 

focusing  of  zero  angular  momentum  ions  vas  used  to  develop  estimates 

of  the  source  brightness  required  to  inject  an  ion  beam  vith  a  given 

current  into  the  transport  channel  in  such  a  vay  that  at  the  focal 

plane,  a  number  density  profile  vith  an  off-axis  peak  located  in  a 

specified  annulus  vould  result.  These  estimates  indicate  that  for 

reasonable  values  of  microdivergence,  the  anode  current  density  that 

is  required  to  deliver  the  desired  amount  of  ion  beam  current  to  the 

2 

target  annulus  is  quite  high  e.g.  -500  kA/cm  for  15  rnrad  and  a  beam 

concentrated  in  the  annulus  6  mm  <  r  <  8  mm.  This  is  for  the  case  of 

an  ion  beam  transported  through  a  3  cm  radius,  11.1  kA  transport 

channel  and  focused  vith  a  300  kA  one  eighth  betatron  vavelength 

focusing  cell.  Hovever,  for  beams  concentrated  in  broader  annuli,  the 

2 

required  current  density  is  substantially  lover  e.g.  -15  kA/cm  at  15 

rnrad  for  r.  =  2  mm  and  r  _  ®  8  mm. 
in  out 

The  effect  of  angular  momentum  has  also  been  considered,  although 
that  analysis  has  not  been  included  here.  One  might  expect  that 
because  angular  momentum  causes  ions  to  stay  avay  from  the  axis,  it 
could  be  helpful  in  producing  the  type  of  annular  focus  considered 
here.  The  detailed  analysis  of  angular  momentum  shovs  that  angular 
momentum  yields  no  positive  benefit,  and  in  fact,  tends  to  make  an 
annular  focus  more  difficult  to  achieve.  The  reason  for  the  lack  of  a 
positive  benefit  may  be  understood  intuitively  by  considering  the 
effect  that  nonzero  angular  momentum  has  on  the  trajectory  that  an  ion 
takes  in  the  radial  phase  space  plane.  Figure  24  shovs  hov  the  radial 
phase  space  trajectory  of  an  ion  is  altered  when  the  ion  is  given 
angular  momentum.  Recall  that  in  the  zero  angular  momentum  case,  the 


phase  space' trajectory  was  approximately  an  ellipse,  but,  as  shown  in 

Fig.  24,  angular  momentum  alters  this  trajectory  by  giving  the  ion  a 

distance  of  minimum  approach  to  the  axis.  The  phase  space  trajectory 

of  the  ion  is  altered  appreciably  from  its  zero  angular  momentum 

trajectory  only  in  the  region  near  the  axis.  This  statement  holds 

true  as  the  ion  passes  through  the  focusing  cell  and  through  the 

ballistic  propagation  region.  Thus,  the  addition  of  angular  momentum 

affects  the  focused  ion  beam  distribution  only  in  the  region  very 

close  to  the  axis.  This  leads  to  the  observation  that  angular 

momentum  does  not  have  an  appreciable  effect  on  enough  of  the  ion  beam 

phase  space  distribution  to  significantly  effect  that  part  of  the 

distribution  that  is  focused  onto  the  annular  target  region.  The  only 

possible  way  in  which  angular  momentum  could  have  a  significant  effect 

2  2 

would  be  if  all  or  nearly  all  of  the  transverse  energy  m(v  +  v0)/2 

r  “ 

available  from  the  ion  diode  were  invested  into  angular  momentum,  and 
this  vould  require  a  diode  that  would  produce  little  or  no  radial 
velocity.  As  such  a  diode  could  not  rely  on  geometric  focusing,  it 
would  have  to  have  a  very  high  source  brightness  to  deliver  a  beam  of 
the  required  current  level. to  the  transport  channel.  Thus,  the 
inclusion  of  angular  momentum  is  not  helpful  for  producing  the  annular 
focused  ion  beam  distribution  considered  here. 
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Appendix 


The  beam  number  density  as  a  function  of  radius  is  obtained  by 
integrating  the  expressions  for  the  beam  phase  space  densi  ty,  Eqs .  (8) 
or  (9),  over  the  velocity  space  coordinates.  For  the  most  part,  this 
amounts  to  doing  integrals  of  the  form 


(Al) 

2  2  2  2 

For  the  case  of  +  ri  <  Yi  +  r2*  corresponding  to  the  phase  space 

distribution  given  by  Eq.  (8),  the  number  density  profile  is: 
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Figure  1:  A  typical  focused  ion  beam  number  density  profile 
exhibiting  the  characteristic  1/r  shape.  The  initial  distribution 
filled  a  rectangular  region  in  the  radial  phase  space  plane  defined 
by  0  <  r  <  2.12  cm  and  0  <  J v^/v^ }  <  0.0598. 


37 


/ALL  STABILIZED.LIGHT* ION  BEAM  TRANSPORT 
AND  FINAL  FOCUSING  SYSTEM 


The  main  components  of  the  ion  beam  production  and 


v  (cm/ns) 


and  focused  ion  beam  radial  phase  space  distribution  for  developing 


the  source  brightness  requirements.  The  intersection  of  the  ellipse 


that  defines  the  ion  beam  phase  space  distribution  vith  the  axis 

must  coincide  vith  the  target  annulus  r.  <  r  <  r 

in  —  out 
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diagram  for  annular  focusing. 
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Figure  8:  The  anode  current  density  for  pinched  focusing  as  a 
function  of  the  focal  length  and  the  microdivergence. 
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Figure  10:  The  minimum  andoe  current  density  for  pinched  focusing  - 
low  current  case. 
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ANODE  RADII  AND  OPTIMAL  FOCAL  LENGTH  (cm) 
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Figure  11:  Inner  and  outer  anode  radii  and  optimal  focal  length  as 
functions  of  the  raicrodivergence  -  low  current  case. 
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ANODE  RADII  AND  OPTIMAL  FOCAL  LENGTH  (cm I 


Figure  12:  Inner  and  outer  anode  radii  and  optimal  focal  lenght  as 
functions  of  the  microdivergence  -  high  current  case. 
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Figure  13:  Source  brightness  as  a  function  of  the  inner  radxus  of  the 
target  annulus  for  pinched  focusing. 
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Figure  14:  The  minimum  anode  current  density  for  annular  focusing 
high  current  case. 
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Figure  17:  Inner  and  outer  anode  radii  and  optimal  focal  lenght  as 
functions  of  the  microdivergence  -  high  current  case. 


ANNULAR  FOCUS 


Figure  18 r  Source  brightness  as  a  function  of  the  inner  radius  of  the 
target  annulus  for  annular  focusing. 


limit  of  the  delevery  efficiency.  The  focused  beam  radius  is  1  cm. 


The  efficiency  is  the  ratio  of  the  arc-length  of  the  darkened  portion 


of  the  quarter  circle  to  the  arc-length  of  the  entire  quarter  circle. 
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Figure  21:  The  delivery  efficiency  as  a  function  of  rin  for  annular 
focusing.  Also. plotted  is  the  delivery  efficiency  in  the  zero 
eraittance  limit. 
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Figure  22;  The  source  brightness  required  to  deliver  the  required  ion 
bean  current  to  the  target  annulus  for  pinched  focusing. 
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